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Abstract 

Non-typical transport phenomena may arise when randomly driven particles 
remain in an active relationship with the environment instead of being passive. 
If we attribute to Brownian particles an ability to induce alterations of the envi- 
ronment on suitable space-time scales, those in turn must influence their further 
movement. In that case a general feedback mechanism needs to be respected. By 
resorting to a specific choice of the particle-bath coupling, an enhanced (super- 
diffusion) or non-dispersive diffusion- type processes are found to exist in generi- 
cally non-equilibrium contexts. 

1 Prerequisites 

A simple Brownian motion, in its canonical (model) manifestations, does not seem 
to hide any surprises. It is our main goal in the present paper to reconsider those 
ingredients of the standard formalism which, while relaxed or slightly modified, would 
lead to conceptually new and possibly exotic (in the sense of being non-typical) features. 

One of possible "defects" of the standard theory is that the Brownian motion is 
incapable to originate spatiotemporal patterns (structures) and rather washes them 
out, even if initially in existence. In this particular context, "active" Brownian parti- 
cles were introduced, 0-0. Most generally, they are supposed to remain in a feedback 
relationship with the environment (heat reservoir, thermostat, random medium, deter- 
ministic driving system or whatever else) through which they propagate. If regarded 
as a thermal bath, the environment is basically out of equilibrium. Obviously, a con- 
crete meaning of being active relies on a specific choice of the model for the thermostat 
(deterministic or random, linear or nonlinear etc.) and the detailed particle-medium 
coupling mechanism. 
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Another "defective" feature of the standard theory is rooted in its possible molecular 
(chaos) foundations where the major simplification is needed while passing from the 
Boltzmann collision theory to the novel kinetic framework set by the Kramers equation, 
cf. ||. In the process, microscopic (molecular) energy- momentum conservation laws 
completely evaporate from the formalism. Then, there is no obvious way how to 
reconcile random exchanges of energy and momentum between the Brownian particle 
and the thermostat with the manifest isothermality requirement. An issue of the 
heat production/removal and resulting thermal inhomogeneities is normally ignored 
(cf. however [|7|-||13| for a discussion of various thermostat models, and the notion of 
thermostatting in non-equilibrium dynamics). 

That was the starting point of our discussion of the origins of isothermal flows and 
the ultimate usage of the third Newton law in the mean to justify the concept of the 
Brownian motion with a recoil (via feedback relationship with the bath), p|. 



In conformity with the metaphor, |TJ|]: "everything depends on everything else", we 



relax the habitual sharp distinction between an entirely passive particle (nonetheless 
performing a continual erratic motion) and its exclusively active (perpetually at rest 
and in equlibrium) driving agency - the thermostat. 

In the present paper we take the following point of view on major signatures of the 
particle activity while moving at random due to external (environmental) reasons: any 
action upon the particle exerted by the environment induces (on suitable space and 
time scales) a compensating reaction in the carrier random medium. In loose terms of 
Refs. [0-0 we attribute to any (even single) Brownian particle an ability to generate 
perturbations of the medium (named a selfconsistent field in Ref. [|lj and perturbations 
of noise in Ref. 0) which in turn influence its further movement. That is a non- 
liner feedback relationship mentioned before: while inducing a random dynamics of a 
particle (action), the medium suffers alterations (reaction), which modify next stages 
of the particle motion (feedback). Such features are obviously alien to the standard 
formalism. 

There is no experimentally reliable way to watch Brownian particles individually 
on time scales below say y^js, nor get any insight into the associated fine detailed 
particle-bath interaction phenomena. The only realistic procedure to quantify observ- 
able features of the Brownian motion is to exploit its familiar Janus face (cf. [l|]). 
Namely, the same mathematical formalism essentially applies to a single particle and 
to a statistical ensemble of identical non-interacting Brownian particles. Then a hy- 
drodynamical analogy can be exploited on both levels of description, f|T|, [R| , once we 
invoke a standard kinetic reasoning (used in passing from the Boltzmann equation to 
gas or liquid dynamics equations). 

Effectively, we need to evaluate (conditional) averages over a statistical ensemble 
of particles, which (in loose terminology that follows Ref. Hl5| , see also |17]]) consti- 
tutes one component (noninteracting Brownian "gas") of a coupled two-component 
continuous system. Another component is the thermostat. 

What are the ultimate kinetic features of the Brownian motion (diffusion process) 
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critically relies on the thermostat model i.e. specific mechanisms of the energy/heat 
exchange due to the coupling between the thermostat and the Brownian "gas". This 
particular issue we shall investigate in some detail in below, by resorting to somewhat 
nontypical (exotic in the Brownian motion context) methods. 

2 Local conservation laws for the Brownian motion 
and Smoluchowski-type diffusion 

It is useful to exploit a standard phase-space argument that is valid, under isothermal 
conditions, for a Markovian diffusion process taking place in (or relative to) a driving 
flow w(x,t) with as yet unspecified dynamics nor concrete physical origin. (In partic- 
ular, such a flow can receive an interpretation of a selfconsistent field generated in the 
environment by Brownian particles themselves, cf. Ref. |lj].) 

We account for an explicit external force (here, acceleration K = F/m) exerted 
upon diffusing particles, while not directly affecting the driving flow itself. Then, 
infinitesimal increments of phase-space random variables read: 

dX(t) = V(t)dt (1) 

dV(t) = p[w(£, t) - V(t)]dt + K(x)dt + (3V2DdW(t) . (2) 

Following the leading idea of the Smoluchowski approximation, we assume that /3 is 
large, and consider the process on time scales significantly exceeding Then, an ap- 
propriate choice of the velocity field w(x, t) may in principle guarantee the convergence 
of the spatial part X(t) of the process to the Ito diffusion process with infinitesimal 
increments : 

dX(t) = b(x, t)dt + V2DdW(t) . (3) 

In this case, the forward drift of the process reads b(x,t) = w(x, t) + ^K(x). Notice 

that the j3~ l K contribution can be safely ignored if we are interested in the dominant 
driving motion. 

Throughout the paper we are interested in Markovian diffusion processes, which 
propagate respectively the phase-space or configuration space probability densities 
(weak solutions of stochastic differential equations are thus involved). In the config- 
uration space variant, we deal with a Markovian stochastic process whose probability 
density p(x, t) evolves according to the standard Fokker-Planck equation 

d t p = DAp - V • (bp) (4) 

and the forward drift is not necessarily of the form dictated by Eqs. (1) and (2). We 
admit here more general forward drift functions, cf. Ref. 0, which do not allow for 
a simple additive decomposition and thus are capable to give account of nonlinearities 
in the particle-bath coupling. 
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One can easily transform the Fokker-Planck equation to the familiar form of the 
continuity equation (hydrodynamic mass conservation law) dtp = — V • (vp) by defining 

v = b — D^£. The current velocity v obeys a local (momentum per unit of mass) con- 
servation law which directly originates from the rules of the Ito calculus for Markovian 
diffusion processes, and from the first moment equation in the diffusion approximation 
of the Kramers theory, |7], [16| : 



d t v + (v- V)v = V(fi-Q) . (5) 

While looking similar to the standard Euler equation appropriate for the lowest 
order hydrodynamical description of gases and liquids, this equation conveys an entirely 
different physical message. 

First of all, for a class of forward drifts that are gradient fields the most general 
admissible form of an auxiliary potential Q(x,t) reads: 

n(x,t)=2D[d t( f>+^(^ + V-b)] . (6) 

Here b(x,t) = 2D'V<f)(x,t). In reverse, by choosing a bounded from below continuous 
function to represent conservative force fields (after taking the gradient) i. e. otherwise 
arbitrary Q, we can always disentangle the above (Riccatti-type) identity with respect 
to the drift field. 

Moreover, instead of the standard pressure term (consider a state equation P ~ 
p a , a > 0), there appears a contribution from more complicated (derivatives !) p- 
dependent potential Q(x,t). It is is given in terms of the so-called osmotic velocity 
field u(x, t) = DV lnp(x, t): 

Q{x,t) = ^u 2 + DV -u . (7) 

An equivalent form of the enthalpy- related potential Q is Q = 2D 2 A J{/ 2 . 

A general expression for the local diffusion current is j = pv = pib-D^). This 
local flow in principle may be experimentally observed for a cloud of suspended particles 
in a liquid. The current j is nonzero in non-equilibrium situations and a non-negligible 
matter transport occurs as a consequence of the Brownian motion, on the ensemble 
average. We thus cannot avoid local heating/cooling phenomena that need to push the 
environment out of equilibrium. That leads to obvious temperature inhomogeneities, 
which are normally disregarded, cf. Ref. 0. 

If the forward drift is interpreted as a gradient of a suitable function (b = 2D'V<f) = 
V$) and we take $(x, 0) as the initial data for the t > evolution), then we have: 

ft = dt$ + i|V$| 2 + M$ . (8) 
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If we decide that the above Hamilton- Jacobi-type equation is to be solved with 
respect to the field t), its solution (and general solvability issue) relies on the choice 
of a bounded from below, continuous function Q(x,t), without any a priori knowledge 
of forward drifts. Viewed that way, Eq. (7) sets limitations on admissible forms of the 
space-time dependence of any conceivable self-consistent field to be generated in the 
bath by Brownian particles (cf. Ref. There is no freedom at all for postulating 

various partial differential equations, if their solutions are to be interpreted as forward 
drifts of Markovian diffusion processes, 



There is also interesting to observe that a gradient field ansatz for the diffusion 
current velocity v = VS allows to transform the momentum conservation law of a 
Markovian diffusion process to the universal Hamilton- Jacobi form: 

n = d t S + ^\VS\ 2 + Q (9) 

where Q(x,t) was defined before as the "pressure/enthalpy" -type function. (That 
form looks deceivingly similar to the standard hydrodynamic conservation law, valid 
for liquids and gases at thermal equilibrium, where instead of Q an enthaply function 
normally appears.) By performing the gradient operation we recover the previous 
hydrodynamical form of the law. 

In the above, the contribution due to Q is a direct consequence of an initial prob- 
ability measure choice for the diffusion process, while Q alone does account for an 
appropriate forward drift of the process, playing at the same time the role of the 
volume force potential (VQ = — contributes to energy-momentum transfer effects 
through the boundaries of any volume). 



3 Moment equations for the free Brownian motion: 
Getting out of conventions 

The derivation of a hierarchy of local conservation laws (moment equations) for the 
Kramers equation can be patterned after the standard procedure for the Boltzmann 



equation, || [16], |TfJ. Those laws do not form a closed system and additional specifica- 
tions (like the familiar thermodynamic equation of state) are needed to that end. In 
case of the isothermal Brownian motion, when considered in the large friction regime 
(e.g. Smoluchowski diffusion approximation), it suffices to supplement the Fokker- 
Planck equation by one more conservation law only to arrive at a closed system. 

To give a deeper insight into what really happens on the way from the phase- 
space theory of the Brownian motion to its approximate configuration-space (Smolu- 
chowski) version, let us consider the familiar Ornstein-Uhlenbeck process (in veloc- 
ity/momentum) in its extended phase-space form. For clarity of discussion, we discus 
random dynamics for one degree of freedom only. 
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In the absence of external forces, the kinetic (Kramers-Fokker-Planck equation) 
reads: 

d t W + uV x W = (3V u {Wu) + qA u W (10) 

where q = D/3 2 . Here (3 is the friction coefficient, D will be identified later with the 
spatial diffusion constant, and provisionally we set D = ksT/mP in conformity with 
the Einstein fluctuation-dissipation identity. 

The joint probability distribution (in fact, density) W(x,u,t) for a freely moving 
Brownian particle which at t = initiates its motion at with an arbitrary inital ve- 
locity Mo can be given in the form of the maximally symmetric displacement probability 
law: 

W(x,u,t) = W(R,S) = [Ax\FG - H*T^.ex P {- GR2 - F g^j^ } (11) 

where R = x- u (l - e"^)/?" 1 , S = u- M e _/3 * while F = f(2/3t - 3 + 4e _ # - e~ m ) 
G = D/3(l - e~ 2l3t ) and H = D(l — e~ m f. 

Marginal probablity densities, in the Smoluchowski regime (we take for granted 
that time scales (3~ l and space scales (D/3 -1 ) 1 / 2 are irrelevant), take familiar forms 
of the Maxwell-Boltzmann w(u,t) = (^tr) 1 ^ 2 ex P(~j^r) an d the diffusion kernel 

w(x,t) = (47r_Dt) -1 / 2 e£]9(— j^t) respectively. 

A direct evaluation of the first and second local moment of the phase-space proba- 
bility density: 

<u>= J du uW{x, u, t) = w(R) [{H/F)R + M e _/3t ] (12) 
<u 2 >= J duu 2 W(x,u,t) = ( FG ~ R2 + ^R 2 )i2nF)-^exp(-^) (13) 



after passing to the diffusion (Smoluchowski) regime, Ul6[ , allows to recover the local 
(configuration space conditioned) moment < u > x = < u > which reads 

<u> x =- = -D v ' 14 
2t w{x,t) 

while for the second local moment < u 2 > T = — < u 2 > we arrive at 

w 

< u 2 > x = (DP - D/2t)+ < u > 2 X . (15) 

By inspection one verifies that the transport (Kramers) equation for W(x,u,t) 
implies local conservation laws: 

d t w + V(<u> x w) = (16) 
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and 

d t {< u > x w) + V x {< u 2 > x w) = -f3 <u> x w . (17) 

By introducing (we strictly follow the moment equations strategy of the traditional 
kinetic theory of gases and liquids) the notion of pressure function P kin (we choose 
another notation to make a difference with the previous notion P of the pressure 
function, cf. VQ = ^f): 

Pkm(x, t) = (<U 2 > x - < U >1)W(X, t) (18) 

we can analyze the local momentum conservation law 

(d t + <u> x V)<u> x =-f3<u > x (19) 

w 

In the Smoluchowski regime the friction term is cancelled away by a counterterm 
coming from -^VPkin so that 

D Vto VP 

(d t + <u> x V)<u > x = - = = -VQ (20) 

2t w w 

where P = D 2 wAlnw. 

For comparison with the notations (local conservation laws) of the previous section 
one needs to replace w(x, t) by p(x, t) and < u > x by v (x, i) in all formulas that pertain 
to the Smoluchowski regime. 

Further exploiting the kinetic lore, we can tell few words about the temperature of 
Brownian particles as opposed to the (posibly equilibrium) temperature of the thermal 
bath. Namely, in view of (we stay in the Smoluchowski regime) Phm ~ {D/3 — %)w 



where D = -^3, we can formally set: 

mmp ' J 

Q = ^^ {kB T-2 l ) < k B T . (21) 
w It 

That quantifies the degree of thermal agitation (temperature) of Brownian parti- 
cles to be less than the thermostat temperature. Heat is continually pumped from the 
thermostat to the Brownian "gas", until asymptotically both temperatures equalize. 
This may be called a " thermalization" of Brownian particles. In the process of ther- 
malization the Brownian "gas" temperature monotonically grows up until the mean 
kinetic energy of particles and that of mean flows asymptotically approach the familiar 
kinetic relationship: 

-(< u 2 > x - < u > 2 x )dx = k B T (22) 

In view of this medium — > particles heat transfer issue, one must be really care- 
ful while associating habitual thermal equilibrium conditions with essentially non- 
equilibrium phenomena, cf. Ref. M for more extended discussion. 
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4 Hydro dynamical reasoning: In search for exotic 



Once local conservation laws were introduced, it seems instructive to comment on 
the essentially hydrodynamical features (compressible fluid/g clS CclSG ) of the problem. 
Specifically, the "pressure" term VQ is here quite annoying from the traditional kinetic 
theory perspective. That is quite apart from the fact that our local conservation laws 
have a conspicuous Euler form appropriate for the standard hydrodynamics of gases 
and liquids. Let us stress that in case of normal liquids the pressure is exerted upon 
any control volume (droplet) by the surrounding fluid. We may interpret that as a 
compression of a droplet. In case of Brownian motion, we deal with a definite decom- 
pression: particles are driven away from areas of higher concentration (probability of 
occurence). Hence, typically the Brownian "pressure" is exerted by the droplet upon 
its surrounding. 

Following the hydrodynamic tradition let us analyze that "pressure" issue in more 
detail. We consider a reference volume (control interval, finite droplet) [ot,0\ in R 1 (or 
A C R 1 ) which at time t G [0,T] comprises a certain fraction of particles (Brownian 
"fluid" constituents). The time rate of particles loss or gain by the volume [a,j3] at 
time t, is equal to the flow outgoing through the boundaries i.e. 

rP 

— d t p(x,t)dx — p(f3,t)v(f3,t) — p(a,t)v(a,t) (23) 



which is a consequence of the continuity equation. 

To analyze the momentum balance, let us allow for an infinitesimal deformation of 
the boundaries of [a,0\ to have entirely compensated the mass (particle) loss or gain: 

[a, /?]->[« + v(a, t) At, (3 + v{(3, t) At] (24) 

Effectively, we pass then to the locally co-moving frame. That implies 

I rP+vpAt r/3 

lim Atl0 — [ p(x,t + At)dx - / p(x,t)dx] = (25) 

At Ja+v a At Ja 

= lim At io——[ p(x, t)dx + At / (d t p)dx + / p(x, t)dx\ — 

At Ja+v a At Ja Jp 

Let us investigate what happens to the local matter flows (pv)(x,t), if we proceed 
in the same way (leading terms only are retained): 

fP+vpAt f P 

/ (pv)(x,t + At)dx- (pv)(x,t)dt~ (26) 

Ja+v a At J a 

rP 

~ -(pv 2 )(a,t)At + (pv 2 )((3,t)At + At / [d t (pv)]dx 
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In view of local conservation laws we have d t (pv) = — V(pv 2 ) + pV(f2 — Q) and the 
rate of change of momentum associated with the control volume [a, (3] is (here per unit 
of mass) 

1 r /3+v At r /3 r/3 

Hm Atl0 — [ (pv)(x,t + At)- / {pv){x,t)] = pV{Q-Q)dx (27) 

lAT Ja+v a At J a J a 

However, VQ = — and P = D 2 pAlnp. Therefore: 

f p V(Q - Q)dx = [ pVttdx - [ VPdx = (28) 

Ja Ja Ja 

= E[Vnf a + P(a,t)-P(P,t) 

Clearly, VQ refers to the Euler-type volume force, while VQ (or more correctly, P) 
refers to the "pressure" effects entirely due to the particle transfer rate through the 
boundaries of the considered volume. 

The missing ingredient of our discussion is the time developement of the kinetic 
energy of the matter flow \{pv 2 ) (per unit of mass) transported through the chosen 
volume. Let us therefore evaluate: 



/ -(pv 2 )(x,t + At)dx- -(pv 2 )(x,t)dt ~ (29) 



~ -^(pv 3 )(a, t)At + ^(pv 3 )(P, t)At + At Jjd t ^(pv 2 )}dx 

We have d t (\pv 2 ) = -\v 2 V(pv) - pvV(Q -Q+ \v 2 ) = -V(\pv 3 ) - pvV(Q - SI). 
Consequently, the time rate of the kinetic energy (of the flow) loss/gain by the volume 
reads: 

Um Am ^-[ f P+V " At hpv 2 )(x,t + At) - f hpv 2 )(x,t)} = [\pv)V(Q - Q)dx (30) 

lAt Ja+v a At I Ja. I Ja 

In the integrand one immediately recognizes an expression (up to the mass param- 
eter) for the power released/absorbed by the control volume at time t. 

In particular, by taking advantage of the identity VQ = ^j- we can rewrite the 

pure "pressure" contribution as follows: — f£ pvVQdx = — J^vVPdx which clearly is 
a direct analog of the standard mechanical expression for the power release (^ = F-v). 

Surely, one cannot interpret the above local averages as an outcome of an innocent 
operation of the random medium upon particles. There is a non-negligible transfer of 
energy and momentum to be accounted for in association with the Brownian motion, 
and a number of problems (when can we disregard temperature gradients ?) pertaining 
to local heating and cooling phenomena suffered by the environment (in terms of local 
averages) should be consistently resolved. 
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5 Implementing the feedback 



In contrast to molecular chaos derivations based on the Boltzmann equation, ||, in 
case of the Brownian motion we have no access into the microscopic fine details of 
particle propagation. Instead, we need to postulate the mathematically reliable form 
of the dynamics, even if ending up with obvious artifacts of the formalism (like e.g. 
non-differentiable paths, or their infinite spatial variation on arbitrarily short time 
intervals, in case of the Wiener process or, even worse, white- noise input). 

Our major hypothesis about an exotic particle-bath coupling pertains to local av- 
erages and is motivated by a fairly intuitive picture of the behaviour of suspended 
particles in thermally inhomogeneous media. Namely, it is well known that hot ar- 
eas support much lower density of suspended particles (can even be free of any dust 
admixtures) than the lower temperature areas. Dynamically we can interpret this phe- 
nomenon as a repulsion of suspended particles by warm areas and an attraction by the 
cool ones, cf. Ref. J7J. 

In the course of our discussion we have spent quite a while on demonstrating that 
a non-trivial energy/heat exchange is completely ignored in the standard approach to 
the Brownian motion (even, if under suitable circumstances one has good reasons to do 
that). Effectively, the model dramatically violates basic conservation laws of physics. 
That derives from the assumption that the thermostat is perpetually in the state of 
rest (in the mean), in thermal equilibrium and thus free of any thermal currents. 

The force/acceleration term (we turn back to the three-dimensional notation) V(fi— 
Q) appears in all formulas that refer to energy /momentum flows supported by local 
averages of the Brownian motion. Their sole reason is the action of the random medium 
upon particles which causes an expansion of the Brownian "swarm" out of areas of 
higher concentration. That however needs a local cooling of the medium (Brownian 
particles are being "thermalized" by increasing their mobility e.g. temperature). On 
sufficiently low time scales that should amount to an instantaneous reaction of the 
medium in terms of thermally induced currents: Brownian particles are driven back 
(attraction !) to the cooler areas i.e. float down the temperature gradients as long as 



they are non- vanishing, [19 



We shall bypass the thermal inhomogeneity issue by resorting to an explicit energy- 
momentum balance information available through local conservation laws for the Brow- 
nian motion. 

If we regard the term V(fi — Q) as a quantification of the sole medium action upon 
particles, then the most likely quantification of the medium reaction should be exactly 
the opposite i.e. V(Q — Q). 

Told otherwise, medium "effort" to release momentum and energy from the Brow- 
nian "droplet" (control volume), at a time rate determined by the functional form 
of V(0 — Q)(x,t), induces a compensating (in view of the heat deficit) energy and 
momentum delivery to that volume needed to remove the thermal inhomogeneity of 
the thermostat. As a consequence, Brownian particles propagate through the medium 
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which is non longer in the state of rest and develops intrinsic (mean, on the ensemble 
average) flows. 



A mathematical encoding of this Brownian recoil principle ||20|| , or third Newton 
law in the mean hypothesis is rather well established. 

The momentum conservation law for the process with a recoil (the reaction term 
replaces the decompressive action term) will read: 

d t v + (v-V)v = V(Q-fl) (31) 

implying that 

d t S + ^\VS\ 2 -Q = -Q (32) 



stands for the corresponding Hamilton- Jacobi equation, cf. plfl , instead of "normal" 
one. A suitable adjustment (re-setting) of the initial data is here necessary, cf. 0. 

In the coarse-grained picture of motion we shall deal with a sequence of repeat- 
able feedback scenarios realized on the Smoluchowski process time scale: the Brown- 
ian "swarm" expansion build-up is accompanied by the parallel counter- flow build-up, 
which in turn modifies the subsequent stage of the Brownian " swarm" migration (being 
interpreted to modify the forward drift of the process) and the corresponding, built-up 
anew counter-flow. 

Perhaps surprisingly, we are still dealing with Markovian diffusion-type processes, 
22| . The link is particularly obvious if we observe that the new Hamilton- Jacobi 



equation can be formally rewritten in the previous form by introducing: 

n r = d t S+^\VS\ 2 + Q (33) 

where Q r = 2Q — Q and Q represents the previously defined potential function char- 
acterizing any Smoluchowski (or more general) diffusion process. It is however Q r 
which would determine forward drifts of the Markovian diffusion process with a recoil. 
Those must come out from the Cameron-Martin-Girsanov identity Q r = 2Q — Q = 

2D[d t <f ) +l(§ 5 + V-b)}. 

After complementing the Hamilton- Jacobi-type equation by the continuity equa- 
tion, we again end up with a closed system of conservation laws. The system is badly 
nonlinear and coupled, but its linearisation can be immediately given in terms of an 
adjoint pair of Schrodinger equations with a potential Q (the imaginary unit i on the 
left-hand-side in below is not an error !), [E2l Ell, |23H. Indeed, 



tdtifj = -DM + ^ (34) 

with a solution represented in the polar form ip = p^^exp^S) and its complex adjoint 
makes the job. The choice of ip{x, 0) sets here a solvable Cauchy problem. Notice 
that, in view of the Schrodinger-type linearization, for time-indepedent Q (conservative 
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forces) the total energy /^(y — Q + Vt)pd z x of the system is a conserved finite quantity. 
We thus reside within the framework of so-called finite energy diffusion processes, whose 
mathematical features received some attention in the literature. In particular, it is 
known that in the absence of volume forces, a superdiffusion appears, while harmonic 
volume forces allow for non-dispersive diffusion- type processes, [[?]]. 
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